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Abstract 

The temperature dependence of the dynamics of water inside microporous activated carbon 
fibers (ACF) is investigated by means of incoherent elastic and quasi-elastic neutron scattering 
techniques. The aim is to evaluate the effect of increasing pore size on the water dynamics in 
these primarily hydrophobic slit-shaped channels. Using two different micropore sizes (~12 and 
18 A, denoted respectively ACF-10 and ACF-20), a clear suppression of the mobility of the water 
molecules is observed as the pore gap or temperature decreases. This suppression is accompanied 
by a systematic dependence of the average translational diffusion coefficient and relaxation 
time (ro) of the restricted water on pore size and temperature. The observed Dr values are tested 
against a proposed scaling law, in which the translational diffusion coefficient Dr of water within 
a porous matrix was found to depend solely on two single parameters, a temperature independent 
translational diffusion coefficient Dc associated with the water bound to the pore walls and the ratio 
6 of this strictly confined water to the total water inside the pore, yielding unique characteristic 
parameters for water transport in these carbon channels across the investigated temperature range. 

Keywords: Supercooled Water, Quasi-Elastic Neutron Scattering, Molecular Diffusion, Porous Carbon 
Fibers. 
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Due to its polar nature, water is considered to be the de-facto solvent of choice for many 
important chemical and biological processes. In the animal kingdom for example, water 
facilitates key essential functions such as the regulation of body temperature, the breaking 
down of nutrients, the activities of vital protein and enzymes, to name a few. These vital 
roles coupled to the extraordinary thermo-physical properties of water, such as its ability to 
expand rather than contract upon cooling, and its known thermodynamics anomalies (heat 
capacity, and isothermal compressibility etc.) continue to drive fundamental research on this 
complex fluid. 

In many common life circumstances however, water is often constrained to some surfaces, 
or restricted within tight cavities or very small voids or cracks. This nano-scopically re¬ 
stricted water undergoes structural perturbations and exhibits reduced molecular mobility 
and thermodynamics properties that are very distinct from those of bulk water [?????? 
]. The slowing down of water dynamics in nano-confinement is rather well documented, and 
has been observed with neutron spectroscopy in a number of porous materials, including the 
hydrophilic silica-based systems such as GelSil [? ], Vycor [? ], SBA-15 [? ], MCM-14 [? ? 
? ? ], and FSM-12 [? ], and the hydrophobic carbon systems such as carbide-derived carbon 
(CDC) [? ], single-wall nanotubes (SWNT) [? ], and double-wall nanotubes (DWNT) [? ? 
]• 

The most commonly investigated pore geometry with neutrons in conhned water research 
is that of the cylindrical type (i.e. silica MCM-41 or carbon SWNT), which is comparatively 
easier to model than non-uniform or other complex shapes. This convenient geometry being 
readily available in various sizes, it facilitates systematic studies and enables direct compar¬ 
ison with proxy models used in molecular dynamics simulations [? ? ]. 

In reality though, natural conhnement of water occurs under various spatial and geo¬ 
metrical restrictions. Understanding the effects of other pore shapes and spacings on the 
dynamics of water is thus of key scientific importance. In the present study, we investigate 
the effect of pore size on the dynamics of water conhned inside carbon channels that are 
primarily slit-shaped. Specihcally, the characteristic diffusive dynamics of water conhned 
in the micropores of two activated carbon hber samples have been investigated at various 
temperatures using neutron spectroscopy. The observed average characteristic relaxations 
(to) and translational dihusion coefficients Dr (associated with the H-sites of water) indicate 
a dihusion that is remarkably slower than in the bulk liquid, with an Arrhenius temperature 
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behavior and characteristic energy barrier Ea that is however similar to the bulk value. The 
diffusion properties can be adequately parametrized using a recently proposed scaling law [? 
] for water transport, leading to two unique water transport characteristics in these slit-like 
carbon pore structures. These observed parameters can be used to predict the translational 
diffusion coefficient of water in these porous carbon samples at any other temperature where 
the law holds, without having to perform additional measurements. 


I. SAMPLE CHARACTERISTICS AND DETAILS 

The KynoF'^* activated carbon fiber (ACT) samples were received from the American 
Technical Trading (ATT), NY. These ACT samples were synthesized from polymeric carbon 
precursors and contain narrow pore size distributions and a large pore volume that is set by 
the degree of activation during synthesis. We use two different ACT samples (ACF-1603-10, 
and ACF-1603-20) for our study. The as-received carbon fibers have average macroscopic 
dimensions of ~ 3 mm length, and ~ 10 /xm diameter. These nominal average dimensions 
were checked against transmission electron microscopy results for consistency (not shown 
here). Previously reported scanning tunneling microscopy (STM) of similar ACF samples [? 
? ? ] has revealed a highly sinuous pore network with predominantly uniform nanometer 
size pore distribution (with a micropore volume V^pore of ~88% of the total pore volume 
Wrotai ~0.4-0.5 cc/g] in ACF-10 [? ? ]). The remainder of the pore volume consists of 
random meso-pores and some ultra-micropores. The relevant sample characteristics are 
listed in Table HI Various other measurements including thermodynamics)? ? ], small- 
angle neutron and X-ray scattering [? ? ] of ACF-10 and ACF-20 have also confirmed 
the interconnected pore structure, made primarily of elongated curvy slit-pores. The pores 
themselves are the voids between curvy but parallel carbon sheets [? ]. The gap between 
two non-flat carbon sheets defines the average pore size. The key point is that the pores 
in ACF are not of cylindrically shaped as in the other well characterized nanotubes, but 
appear to made of irregular slits pores, as indicated by the above cited STM work. This 
different pore structure offers a new platform for investigating fluids in confined geometries 
in non-cylindrical pores. For the current ACF-10 and ACF-20 samples, the average pore 
size is estimated to be respectively (12T0.5) and (18±0.6) A, based of previously reported 
N 2 adsorption isotherm measurements [? ? ] using the Dubnin-Astakhov equation [? ]. 
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TABLE I: Characteristic size and volume of the micropores in ACF, corresponding BET surface 
area (from Ref. [? ]) and hydration level h (weight%), as measured from the relative weight 
change of ‘dry’ sample after 24 hours exposure to a humid atmosphere. E^pore is the open volume 
of micropores, and Vrotai the total available pore volume in the sample. 


Sample 

p-pore size (A) 

Vfipore 

Vrotai ^ 

Surface Area (m^/g) 

h (weight%) 

ACF-10 

11.97 

88 

785 

21 

ACF-20 

17.69 

75 

2247 

23 


Since the average pore sizes are all below 2 nm, we refer to the samples as microporous 
materials. The nominal sample specihc surface areas, tabulated in Table [H were conhrmed 
at the time of the synthesis by the manufacturer by iodine number testing, and correlated 
well with BET calculations in this range [? ]. Iodine is a common standard adsorbate used 
in industry to estimate the adsorption capacity of carbon samples. 

To prepare for the neutron experiments, we outgassed the as-provided ACF samples for 
36 hours at 473 K in a vacuum oven to remove all of the bulk-like water, and most of 
the surface water that was originally present in the as-received sample. We then exposed 
respectively about 2 g of each so-dried sample to a humid atmosphere in a desiccator for 
several hours. The hydration level reached in each case after about 24 hours exposure is 
indicated in Table HI These values were based on the relative weight change of each sample. 
These somewhat important hydration levels signal the presence of hydrophilic groups (such 
as oxygenated sites for example) in otherwise totally hydrophobic samples. Since the quoted 
hydration amounts are relative to the ‘drying’ conditions set above, it is important to note 
that the diffusive dynamics reported in this work are those of all water molecules present 
inside the porous carbon network. The neutron being primarily sensitive to hydrogen atoms, 
the present measurements yield the characteristic relaxations of all conhned water molecules. 
The hydrated samples were subsequently each loaded unto two concentric A1 cylinders (with 
a 2 mm gap between them) to minimize multiple scattering. The subsequently indium-sealed 
containers were anchored to the copper huger of a close-cycle refrigerator (CCR) stick which 
allowed to control the sample temperature for the present measurements between 50 and 
300 K. 
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FIG. 1: (Color online) Observed quasi-elastic neutron scattering intensity I{Q,E) on the BASIS 
instrument, as a function of energy transfer E for water in ACF-20 (average slit size d ~ ISA) at 
selected momentum transfer Q and temperature T = 250 K. Although the actual energy window 
investigated is symmetric in energy (±120 /reV), only the neutron energy gain side is shown for 
display purposes. The relative errorbars are smaller than the symbols. 

II. NEUTRON SCATTERING MEASUREMENTS 

In this section, the technical details concerning the neutron measurements are presented. 
The neutron data were all collected on the backscattering spectrometer (BASIS) at Oak 
Ridge National Laboratory (ORNL), USA [? ]. This unique indirect geometry neutron spec¬ 
trometer has an excellent energy resolution T^ = 1.75 fieV (Half-Width-at-Half-Maximum 
or HWHM) at the elastic line, and covers a momentum Q and energy E transfer range, 
respectively 0.3 < Q < 2 A“^, and —120 < E < 120 /reV, enabling access to time and 
length scales in the range I =3-22 A and t from ~5 up to 1000 ps. To show the quality of 
the data collected on BASIS, representative S{Q,E) spectra of water adsorbed in ACF-20 
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FIG. 2: (Color online) Temperature and pore size dependence of the energy-integrated elastic 
intensity I{T) = dE St{Q,E) of water confined in carbon fibers, at a selected momentum 
transfer Q of 1.3 A, where F^ is the instrument energy resolution width (HWHM=1.75/reV). The 
onset of observable diffusive dynamics on BASIS is marked by the departure from a monotonically 
varying slope in I{T) at low temperatures to a more rapidly changing slope at the higher tempera¬ 
tures. The solid lines are model fits to data, corresponding to an activation energy of Ea =16 and 
21 kJ/mol, for ACF-10 and ACF-20 respectively. The errorbars are smaller than the symbols. 

at temperature 250 K is shown in Fig. [H 

A. Incoherent Elastic Response 

To gain preliminary insights into the molecular dynamics of water in ACF and determine 
the appropriate temperature range within which the diffusive dynamics become observable 
on BASIS, it is useful to investigate the temperature dependence of the elastic peak. This 
is achieved by integrating the peak intensity of rapidly collected runs (10 mins per point) of 
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S{Q, E) around E = Q (—< E < F^) at each temperature from say 100 to 300 K, in steps 
of 5-10 K depending on the temperature region. Fig. [2]shows the result of such an analysis for 
both ACF samples, at a selected Q of 1.3 A. The resulting energy-integrated elastic intensity 
/(T) is normalized to the lowest temperature data taken at Tq = 50 K so that /(T)//(To) 
equals unity at T = Tq. The I{T) signal decreases with increasing temperature, indicating an 
increase in mobility of the H-sites of water, and consequently of water itself. This behavior 
can be easily understood on the basis of an I{T) that is modulated by a Debye-Waller 
coefficient; i.e. proportional to exp [—(5 ^(m^(T))/ 3] [? ? ] at each Q, where {v?{T)) is the 
mean square displacement (MSD) associated with the hydrogen atoms in water. At low T, 
the MSD arises primarily from harmonic vibrations, increasing monotonically with increasing 
T. At high T where a n h armonic vibrations become important, its accurate determination 
requires a subtraction of the harmonic phonon contributions (which can be estimated from 
the MSD at low T by extrapolation). 

In any case, the total elastic area I{T) integrated over the instrument resolution window 
would effectively increase as T is reduced, converging to some constant value when the 
diffusive dynamics become resolution limited, which is observed here to be around 150-180 
K. Similarly, when the dynamics become too broad to be observed at the high temperatures 
(above ~ 290 K here), I(T) flatten to background levels. These levels are dominated by 
higher energy dynamical processes that fall outside the instrument dynamics range. The 
intermediate region between these two limits where the intensity drops relatively quickly is 
associated with the relevant molecular diffusion. Water crystallization on the other hand, 
if present, would manifest itself along the I{T) curve as a sharp drop (similar to a step- 
like hrst order transition), which is not observed here at any Q (i.e. at any length scale 
probed). The scattering from solid ice is purely elastic. Our present observation suggests 
that water molecules in the microporous ACF remain mobile well below bulk water freezing 
temperature. 

To get a preliminary estimate of the energies associated with the thermally activated water 
dynamics suggested by the drop in I(T) at various length scales, the diffusive component 
of S{Q,E) can be modeled with a single Lorentzian function Lt{E) = to a hrst 

approximation. Assuming an Arrhenius temperature dependent Tt = Fqo exp(—£'^/i?T)), 
the energy-integrated I{T) for < F^ can be explicitly evaluated at each Q as follows, [? 
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where is the HWHM of the instrument energy resolution (~1.75 /ieV), Too that of Lt{E) 
when T —)■ cxo, and Aq a fraction of the observable molecules that do not participate to 
the diffusive process. Eq. [T]has thus three adjustable parameters: Aq, Too, and Ea_. The 
solid lines in Fig. 0] represent the hts obtained at Q =1.3 using Eq. [H These hts 
yield average Ea values between 13 and 16 kJ/mol for ACF-10 and 16 and 21 kJ/mol for 
ACF-20 for the Q-range probed (0.5 < Q < 1.3 A“^). At the higher Q (> 1.5 A“^), the data 
analysis is complicated by the coupling between translational and rotational modes. These 
Q values are thus not included in the present analysis. Estimate of the energy barriers Ea 
associated with the long range translational diffusion can be best inferred from the lowest 
Q = 0.5 A“^ data (spatial scale up to 21 :jQ ~ 13 A). Similarly, Ea associated with localized 
dynamics (sphere of radius ~5 A) is provided by the largest investigated Q of 1.3 A“^. The 
observed Ea within this small Q interval increases with increasing Q, as indicated in the 
inset of Fig. |2l There is a relatively broad distribution of Ea values, whose observable 
limits are more less set by the lowest and highest accessible Q on the neutron spectrometer. 
This distribution of Ea, as summarized in Table UTl suggests that a single Lorentzian alone 
cannot fully describe the entire quasi-elastic neutron (QENS) spectra. To map-out the entire 
observable Ea values, I{T) should be analyzed at each accessible Q value. Eq. [T] provides 
only a rough estimate of Ea at the corresponding length scale probed. A more sensitive 
and accurate approach for determining Ea for the observed diffusive process is to £t the 
neutron spectra with an appropriate model at each wavevector independently (in either the 
energy or time domains), and investigate the temperature dependence of the corresponding 
translational diffusion coefficient or relaxation times, as we have done below. The parameter 
Afj decreases with increasing Q, varying for example from 0.46 ai Q = 0.5 A“^ to 0.27 at 
Q= 1.3 A-^ in ACF-10. 
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TABLE II: Estimate of the activation energy Ea (in units of kJ/mol) of water in ACE, as obtained 
from fits of Eq. [T]to the integrated elastic intensity I{T). Within the Q-range investigated (0.5 
< Q < 1.3 A“^), these estimates provide lower and upper limits for Ea as indicated in the inset 
of Fig. [2j The low Q yields the characteristic Ea value for long range translational diffusion (up 
to ‘I'KjQ — 13 A) and the high Q gives the characteristic value for the localized dynamics (at a 
length scale of ~5 A). The rather broad distribution of Ea values indicates the limitation of the 
method used. To more accurately determine Ea, an alternative approach based on the temperature 
dependence of the diffusion coefficient (as done below) is needed. 


Q (A-Q 

ACF-10 

ACF-20 

0.5 

13.6 

16.3 

1.3 

16.2 

21.2 


B. Quasi-Elastic Neutron Scattering (QENS) 


Having determined a snitable temperatnre range for investigating the diffnsion of water 
confined in the hydrated ACF samples on the spectrometer, high statistical qnality quasi¬ 
elastic neutron (QENS) data were collected at three temperatures: 280, 250, and 230 K. 
The relevant momentum transfer range for this comparative QENS study is limited to 0.5 < 
Q <1.3 A“^, AQ = 0.2 A“^, to avoid the undesirable influence of the coherent contributions 
of carbon at low and high Q’s, and of the faster rotational motions of water that complicate 
the interpretation of the data. As indicated above. Fig. [1] shows representative spectra as a 
function of Q for water in ACF-20 at 250 K. The hgure also depicts the resolution function 
at Q = 1.3 A“^ taken with the same sample cooled down to 50 K (dashed black line). 
Each QENS spectra was normalized against the same vanadium run to correct for detector 
efficiency, and subsequently Fourier transformed to a self-intermediate scattering function 
I{Q,t) using. 




dE exp{iEt/h) S{Q,E) 


( 2 ) 


— Err 


where Em is ~120 /reV. The instrument contributions were removed by dividing the resulting 
I{Q, t) at each temperature by that obtained at 50 K data. The intrinsic Iin{Q, t) were then 
all htted in the time domain over a limited but reliable t range, 30-750 ps. In principle. 
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the accessible times on BASIS is a low as 5 ps and as high as 1000 ps, given the dynamics 
range and the energy resolution of the instrument. The narrower time domain over which 
the present analysis was done is one for which the translational dynamics of water are 
conveniently relevant on BASIS, and (2) where known systematic errors associated with 
complex Fourier transform of neutron data have minimal effects. Further details regarding 
the Fourier transform methods used here can be found in Ref. [? ]. 


III. DATA ANALYSIS 


Fig. [3] compares the resolution deconvoluted intrinsic Iin{Q,t) for Q = 0.9 at two 
selected temperatures. These resolution independent Iin{Q, t) could be adequately described 
by a Kohlrausch-Williams-Watts (KWW) stretched exponential model [? ] plus a time- 
independent elastic component, 


liniQ, t) = Ao (1 - Ao) exp 



( 3 ) 


where Aq denotes the fraction of immobile water molecules (most notoriously known as 
elastic incoherent structure factor or EISF), r is the relaxation time, and (3 the stretching 
exponent. Our Aq parameter is offset from the true EISF because ‘dry’ sample contribution 
has not been subtracted here. In our recent investigation of water in ACF-10 [? ], we 
found its dynamics to be more reliably captured by the KWW model, rather than the two 
exponential model. This KWW model is suggestive of and consistent with a heterogeneous 
diffusive dynamics. Typical KWW hts can be seen as black solid lines in Fig. |3l As can be 
appreciated, the hts clearly capture the experimental data quite reliably over a wide time 
range, from about 30 ps to about 750 ps. From these hts, the temperature and wavevector 
dependence of the three adjustable parameters Aq, r and (3 are determined. By analyzing the 
Q-dependence of these parameters at each temperature, the dihusion coefficient, and average 
relaxation time (yg) for the conhned water can be determined. Assuming a distribution of 
relaxation times at temperature T, the mean value of (rg) at that temperature is thus the 
area dt = (r//3)F(l//3), where F(x) is the gamma function. 
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FIG. 3: (Color online) Intrinsic Intermediate scattering function Iin{Q,t) of water confined in 
activated carbon fibers, as obtained from Fourier transforming the observed S{Q,E) data, and 
removing resolution contributions at a selected Q of 0.9 A, as explained in the text. The top panel 
compares the results derived for water in ACF-10 (open circles) and ACF-20 (closed circles) at 
temperature T = 280 K. The bottom panel makes a similar comparison of the data collected at 
230 K. The solid black lines are model fit using the stretched exponential model, described in the 
text. The errorbars associated with Iin{Q,t) are less than 1% within the time range shown. 
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TABLE III: Confining cage size a*, and fraction /j of immobile of hydrogen atoms, obtained from 
fits of Eq. S]to the elastic incoherent structure factor shown in Fig. [H 


r(K) 

aio (A) 

^20 (A) 

fio (%) 

f20 (%) 

280 

4.8(2) 

4.7(1) 

47 

29 

250 

4.4(1) 

4.4(3) 

43 

25 

230 

4.0(1) 

4.6(1) 

42 

25 


IV. RESULTS AND DISCUSSION 

The Q-dependence of the Ao parameter (from Eq. [3]) associated with each ACF sample 
is displayed in Fig. IHfor the three investigated temperatures. This parameter introduced 
above is a measure of the fraction of H atoms that do not contribute to the diffusion process 
within a cage of certain dimension (~ 271/Q). With this in mind, it is evident (based on 
Fig. H]) that (1) the fraction of mobile water molecules is largely temperature independent 
(except perhaps in ACF-10), and (2) that at short distances (high Q), approximately 20% 
more water molecules will be contributing to the observed dynamics in the large pores of 
ACF-20 than in the narrow ACF-10 pores, in agreement with Fig. [2l To estimate the size 
Oj (where i =10, 20 for ACF-10, and ACF-20 respectively) of the conhning cage from Aq, a 
model fit is necessary. Using the following generic expression EISF for spherically confined 
motion, [? ], 

-4o = /i + (1 - /.) W 

where ji is the spherical Bessel function, the parameter a* and the actual fraction /* of 
immobile water molecules outside the cage are extracted at each temperature, as summarized 
in Table Hm for both samples {i = ACF-10, ACF-20). These values confirm the larger fraction 
of immobile water molecules in ACF-10 inferred from Fig. 01 and reveal an average conhning 
cage radius for water in ACF-20 of about 4.5 A that is unaffected by temperature (less than 
2% change on cooling from 280 K to 230 K). In contrast, the corresponding cage size in 
ACF-10 shrinks by as much as 16 % for the same temperature change. 

Eq. |3] yields a stretching exponent (3, as well as a relaxation time r that are both 
temperature and wavevector dependent. While the inhuence of these variables on /3 is 
rather marginal compared to that on r and generally kept hxed to some average value in 
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FIG. 4; (Color online) Temperature and momentum transfer dependence of the observed (offset 
elastic incoherent structure factor) for ACF-10 (open circles) and ACF-20 (solid squares). 

KWW fits, we here chose to let it vary with r, and Aq. We subsequently computed the 
average relaxation time (rg). The variation of (3 with Q for water conhned in the ACF 
samples at 250 K is illustrated in Fig. [5l Within the temperature and Q-range investigated, 
the observed (3 values fall in the range of 0.6 < (3 < 0.8, in excellent agreement with previous 
hndings.[? ? ? ] 

From the Q-dependence of (rg), it is possible to determine the nature of the diffusion 
process, whether for example it is translational (quadratic in Q) or rotational [Q invariant) in 
character. The computed values for all the three investigated temperatures are displayed in 
the form of 1/(rg) versus in Fig. [6l The resulting Q-dependence suggests a translational 
jump diffusion, with a quadratic behavior at low Q’s, and a saturation to a hxed value at 
high enough Q. This model is given by, 

^ ^ DrQ^ 

(r^) l + DrQ^{T,) 
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FIG. 5: (Color online) Best fit stretching exponent parameter /3 as a function of Q at T = 2hQK 
for ACF-10 (open circles) and ACF-20 (solid squares), determined within 0.5% precision. The 
solid and dashed lines represent the Q-averaged (/3) values for ACF-10 (0.68) and ACF-20 (0.71), 
respectively. In the present work, [3 for each Q and T were used instead of (/3), as is often the case. 

where Dr corresponds to the average translational diffnsion coefficient of water inside the 
pores (associated with the H-sites). The lines in Fig. [6] represent the best hts of Eq. [5] to 
the experimental data, from which the temperatnre dependence of (tq) and Dr are extracted 
and examined below. 

Fig. [ 7 ] compares the temperatnre dependence of (tq) for water in ACF-10 and ACF-20 
with that of bnlk water, plotted on a logarithmic scale. Within the instrnmental precision, 
onr observed relaxation time of water (tq) changes by a factor of ~6 between 230 K and 280 
K in ACF-10, compared to only 4 in ACF-20, as indicated by the different slopes of the lines 
in Fig. [71 Since shorter times are associated with faster dynamics, one can immediately 
conclnde that while the dynamics of water becomes progressively slower as the conhning pore 
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FIG. 6: (Color online) Inverse of the average relaxation time as a function of and temperature. 
The solid squares are the characteristic values obtained for ACF-20 and the open circles those for 
ACF-10. The Colored lines are the best representative fits using the jump diffusion model discussed 
in the text. The associated errorbars are between 2-9%. 

size is reduced, the corresponding (tq) values between jump-sites, are consistently larger in 
the smaller pores at all temperatures. 

A subsequent investigation of the translational diffusion coefficient Dr shown in Fig. [8] 
conhrms the slowing down of the dynamics as the conhning pore dimension is reduced. 
The behavior of Dr indicates a progressively decreasing mobility of water molecules as the 
conhning space becomes tighter. At 280 K for example. Dr is respectively 40 and 30% lower 
in ACF-10 and ACF-20 than in the bulk liquid at the same temperature. As the temperature 
drops. Dr departs signihcantly from its bulk value, decreasing by as much as 77% in ACF-10 
and 70% in ACF-20 at 230 K. 

Both Dr{T) and {tq){T) can be ht to an Arrhenius exponential decay ~ Aexp(±i7^/i7T), 
yielding the characteristic activation energy Ea associated with each. The resulting hts are 
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TABLE IV: Activation energy E^i (kJ/mol) inferred from various parameters associated with 
water diffusion, where i =10, 20 for water in ACF-10 and ACF-20 respectively. Corresponding 


bulk water Ea = 

value derived from the data in Ref. [? ? 

is shown for comparison. 

Source 

(to) 

Dr (10“®cm^/s) 

Ea^o (kJ/mol) 

15.5 

21.3 

Ea 2 o (kJ/mol) 

24.1 

22.2 

(kJ/mol) 


22.4 


indicated by the dashed lines in Figs. [7] and [HI with Ea values that are summarized in 
Table HVl The Ea values inferred from Dr{T) are approximately 21-22 kJ/mol, comparing 
well with previously reported values in porous media [? ? ]. Contrary to those deduced 
from (ro)(T), they also appear to be more uniform across the two samples, and comparable 
to those of bulk water [? ? ]. Although this estimated bulk water Ea (~22.4 kJ/mol) is 
signihcantly larger than that reported in Ref. ?? for strictly rotational relaxations (E^ =7.7 
kJ/mol), we hnd it to agree rather well with recent Endings in which the corresponding 
molecular motions were not considered to be strictly rotational, with Ea =22-24 kJ/mol [? 
]. Compelling arguments of why the localized water dynamics cannot be considered to be 
purely rotational in character can be found in Ref. [? ]. 

Because the Ea derived from the temperature dependence of (tq) in Fig. [7|do not exhibit 
a clear systematic pore size dependence, we focus instead on the Ea that dictates the 
temperature dependence of D^iT) (see Table HVl) . The results suggest that while the overall 
diEusive dynamics of the water molecules slow down as the pores size or the temperature 
is decreased, the activation energy for H-bond breaking in water confined in porous ACF 
remains globally unchanged from its bulk value. 

In a recent molecular dynamics simulation (MD) work, Chiavazzo et al. [? ] reported a 
‘universal’ scaling law that can be used to interpret water transport in confined geometries. 
This scaling depends primarily on a single parameter, 6 which is the ratio between the 
strictly confined water (or water influenced by the pore walls) to that of the total water 
in the confining pore. Obviously the larger the pore, the smaller this parameter is. They 
tested and validated their predictive model against some 60 cases, ranging from water in 
nanoporous silica to hydration layers in proteins. Based on this proposed scaling model, the 


16 










FIG. 7: (Color online) Mean relaxation time (tq) of water confined in ACF-10 (open circles) and 
ACF-20 (solid circles) plotted as a function of 1000/T. Data for bulk water [? ? ], observed 
at comparable temperatures are also shown. Dashed lines are fits to an Arrhenius temperature 
dependence. 

translational diffnsion coefficient Dj. of confined water can be written at any temperatnre 
as, 

Dr.{T)=eD,+ {l-e)D,^ik{T) (6) 

where Di,uik{T) is the bnlk water translational diffnsion at temperatnre T, and Dc and 6 
the temperature insensitive parameters, respectively the translational diffusion coefficient 
associated with the strictly conhned water (i.e. the water that is most strictly affected by 
conhnement), and 6 the ratio between this water and the total water in the conhning media. 
In practice, Di,uik{T) can be obtained from the literature, but and 6 would be sample 
specihc and geometry dependent. 

Using Eq. |6l it is thus possible to determine the characteristic transport parameters Dc 
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TABLE V: Observed temperature independent characteristic parameters of water confined in ac¬ 


tivated carbons; the fraction 9 of water molecules closest to the pore walls ( ^°otat )’ corre¬ 

sponding diffusion coefficient Dc 


Parameters 

P — ^hnii.n.d 

Dc (10 ®cm^/s) 

ACF-10 

0.41 

0.35 

ACF-20 

0.28 

0.14 



Q 



FIG. 8: (Color online) Net translational diffusion coefficient as a function of lOOO/T for ACF-10 
(open circles) and ACF-20 (solid squares), as compared to that of bulk water (solid diamonds). 
Dashed lines are fits to an Arrhenius behavior ~ 

and 6 for water in the present ACF samples. Fig. [9] shows the observed D^iT) valnes in 
ACF versns known bnlk water data (DhuikiT)). The dashed lines are hts of Eq. |6]to the 
experimentally observed values. From the slopes and intercepts of these lines, 6 and Df. 
can be extracted. The observed values for each sample are listed in Table |Vl Fig. [9] also 
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shows the limiting situations by the strictly conhned water (6* = 1) and pure bulk water 
(6^ = 0). The grey shaded area is inaccessible and corresponds to the unphysical condition 
DriT) > DhuikiT). 

In our opinion, the temperature region over which Eq. [H]is applicable may depend on both 
the sample and the characteristics of the neutron instrument used (resolution and dynamics 
range). Based on the current data, we estimate this range to be between T^m ~ 220 K 
and Tjjiax 290 K on the BASIS spectrometer. The maximum value is primarily set by the 
accessible time window on the instrument, and the lower limit varies based on a combination 
of instrumental resolution, and the physics of water itself (conhning pore size, possible phase 
transition related to thermodynamics anomalies [? ] etc.). It may however be more prudent 
to set Tmin above at least 250 K, the lowest temperature at which bulk supercooled water can 
be easily achieved and for which neutron scattering work has been fairly well documented 
[???]. The law requires thus that the temperature range considered be one in which 
Dc and 6 do not vary with temperature. This suggests a temperature interval in which the 
effective potential well in which water sits in the conhning pore remains greater than the 
thermal kinetic energy. 

In the current study, the translational diffusion coefficient for bulk water at 230 K shown 
in Fig. |9] (lowest value where Dr for ACF-10 and ACF-20 appears to overlap) was estimated 
by extrapolating the temperature dependence of D}yuik{T) above 250 K down to 230 K. This 
approximation is off course by no means exact since Dr of water is known to be rather super- 
Arrhenius [? ], but it is one that yields a reasonable estimate for the present relatively small 
temperature interval. Excluding the data at 230 K does not affect our Endings. In fact, the 
remarkable observed linearity in Fig. |9] suggests that a single temperature point would be 
sufficient to predict the overall temperature behavior, and the characteristic parameters Dc 
and 9. This observation warrants further investigation. 


V. CONCLUSIONS 

In the present study, the molecular dynamics of water adsorbed in microporous activated 
carbon hbers have been investigated for two different pore sizes, and compared with bulk 
water dynamics. The quasi-elastic neutron scattering data of ACF-10 ( 12 A) and ACF-20 
(~ 18 A) reveal a retardation of the water dynamics when either the pore dimension or the 
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230 K 250 K 280 K 



FIG. 9: (Color online) Translational diffusion coefficient {Dr) of water restricted within the pores 
of activated carbons at the three temperatures investigated plotted against the corresponding 
values for bulk water {Dhuik)- The lowest value (230 K) is based on an extrapolation of 

the Arrhenius behavior of shown in FiglSl The dashed lines represents the scaling behavior 

Dr{T) = ODc + (1 — 0)Dbuik{T) proposed in Ref. [? ], where Dc is a temperature independent 
diffusion coefficient associated with the strictly confined water (i.e. most inffuenced by the pore 
walls), and 9 the ratio between this water and the total water inside the pores. The grey shaded 
area is not physically accessible since Dr{T) cannot exceed the bulk value. The characteristic 9 
and Dc parameters for water diffusing in ACF-10 and ACF-20 are summarized in Table IVl 
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temperature is reduced. The observed translational diffusion coefficients can be adequately 
described by a recently proposed scaling law for water diffusion in nano-confined geometries, 
yielding intrinsic characteristic parameters capable of predicting water mobility in each of 
these two specific porous carbon materials at any temperatnre over a wide temperature range 
where the law remains valid. The significance of this scaling law to neutron experimenters, 
who far too often struggle to complete a series of temperature scans within a limited beam¬ 
time allocated by the neutron facilities, cannot be underestimated. The present findings 
snggest that for a given nanoporous media, knowledge of the water dynamics at 2 temper- 
atnres (3 to confirm accuracy, as we have done here) is all that is required for predicting 
its diffnsion at other temperatnres. To farther test the universality of this law, additional 
work is however needed. Fnture research conld involve other well studied porous media for 
which the dynamical properties of water are well known (MCM-41 for example). It will 
also be interesting to farther investigate the temperatnre range over which the method can 
be reliably applied using other neutron spectrometers with different energy resolution and 
dynamics range. 
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